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1.  INTRODUCTION 

Christaller  (1933)  postulated  a regular  hexagonal  lattice 
as  his  model  for  the  distribution  of  towns  over  a homogeneous 
area.  This  is  known  as  Central  place  Theory.  This  is  con- 
trasted  against  the' hypothesis  that  the  spatial  pattern  of 
towns  in  the  area  is  random.  A survey  on  this  and  related 
topics  is  given  by  Cliff  & Ord  (1975)  as  well  as  an  extensive 
bibliography.  See  also  Berry  & Pred  (1961)  and  Ripley  (1977). 

In  an  ideal  situation  Central  Place  Theory  would  predict  a 
pattern  as  given  in  Fig.  1 for  11  central  places.  As  shown, 
the  region  can  be  covered  by  hexagons  of  the  same  size.  We 
will  require  the  following  definitions. 
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Spatial  pattern  under  Central  Place  Theory 
with  11  towns  (•  denotes  a town) 


( i)  Dirichlet  Cells  (Thiessen  Polygons) 

2 

Let  X,,  be  the  co-ordinates  in  R of  n 

~i  ~2  ~n 

central  places.  The  ith  Dirichlet  cell,  n(x^),  is  defined 
such  that  for  all  x C 11  (x^), 

1 1 x-x^ll  < 1 1 x-x^  II  , j / i . 

The  boundary  of  the  cell  defines  the  corresponding  Thiessen 
Polygon  and  these  cells  or  polygons  form  a non-overlapping 
covering  of  the  plane. 

( ii)  Delaunay  Triangles 


These  are  formed  by  connecting  those  points  whose  Theissen 
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poiygons  have  a common  side.  Provided  that  4 or  more  DirichJol 
cells  do  not  meet  at  any  one  point,  the  Delaunay  triangul.ition  will 
also  form  a non -over lapping  covering  of  the  plane.  For  4 cells 
to  have  a common  vertex,  the  4 centres  would  have  to  lie  at  the 
vertices  of  a square  or  rectangle  and  since  we  are  in  continuous 
space,  this  probability  can  be  neglected. 

As  a numerical  example,  we  examine  the  data  of  Gaile  & Burt 


(1977,  -tnpublished)  . They  consider  a map  of  44  Central  Places  in 
6 counties  in  Iowa,  namely  Unwin,  Ringgold,  Clarke,  Decatur,  Lucas 
and  Wayne  Counties.  Thiessen  polygons  were  computed  and  the  interior 


angles,  a^.  of  the  63  resulting  Delaunay  triangles  were 

obtained.  Of  course  a^+  a^+  a^=  180°.  Define  *^(2)'  *^(3) 


so  that 


Table 


gives 
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2.  RANDOMNESS  HYPOTHESIS  AND  MILES'  DI^iNS  ITY 
2.1  Miles ‘ Result 

When  Central  Place  Theory  does  not  hold,  the  T)i lesson  polyqons 
become  irregular  and  the  Delaunay  triangles  follow  a more  random 
pattern.  The  one  procedure  is  then  to  generate  the  points  in  the 
plane  fay  the  homogeneous  Poisson  process.  This  process  has  been 
studied  in  Miles  (1970)  with  the  following  results. 

Let  X be  a 'nice'  domain  in  Rf  e.g.  the  interior  of  a circle. 

Let  q > 0 be  a real  number  and  let  X(q)  be  the  circle  of  radius  q. 

Consider  a Poisson  process  P in  R^  of  constant  parameter  p, 

and  for  those  points  lying  in  X(q)  consider  the  associated  Delaunay 

triangulation.  A triangle  with  vertices  x^,  x^  can  be  re- 

parameterised  in  terms  of  x^^,  R,  6^^,  6^,  where  R is  the 

circumradius  and  6^  is  the  angle  formed  between  the  horizontal 

through  the  centre  of  the  circumcircle,  and  the  line  joining  x^ 

to  the  centre. 

Therefore,  in  the  notation  of  Section  1, 

sina^  = lsinis(e2-  ^ |sin^5(e^-  G^)!.  (0,ti). 

Let  R and  (aj^,a2)  be  asymptotically  independent  with 

g(R)  = 2(np)^R^  exp  (-npR^),  R > 0,  p > 0,  (2 

8 

sin  sin  sin  (cl^+  a^) , a^>  0,  (i^>  0,  a^+  a^<  n,  (2 
and  let  B R^  be  an  open  set  in  (R,  aj^,a2)  space.  Miles'  result 


can  be  stated  as 
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Number  of  dA's  in  X(q)  such  that  (F<,  > d 

Number  of  DA's  in  X(q) 

a.s. 

-*  J g(R)  f (a^,  a2)dRda^da2  (2.M 

B 

where  DA  stands  for  a Delaunay  triangle. 

The  convergence  is  almost  sure  with  respect  to  the  Poisson  process 
P as  q . In  practice  we  have  a set  of  data  with  say  N^  points 

in  a region  X,  and  we  would  like  to  know  if  those  points  could  be 
the  realisation  of  a Poisson  process  and  hence  if  Miles'  results  are 
j applicable  to  the  finite  situation.  It  seems  plausible  that  a process 

closer  to  the  data  would  be  obtained,  if  a priori,  we  condition  on 
there  being  N = points  in  the  region.  Of  course  now  the  process 

is  no  longer  Poisson  and  it  is  natural  to  ask  if  the  asymptotics  of 
Miles  are  valid  for  this  conditioned  process. 

For  N 2 N^,  define  a conditioned  Poisson  process  in  X(q)  of 
N points.  AS  N varies  these  processes  are  defined  on  different 
probability  spaces,  so  it  is  not  possible  to  state  "almost  sure'  re- 
sults, however,  we  still  have  convergence  in  distribution  to  (2.3). 


2 . 2 Verification 

Miles'  result  says  nothing  about  the  distribution  of  ( i ^ 
i = l,...,n.  We  hypothesize  that  the  independently 

distributed  as  f(aj^,a2)  given  by  (2.2)  for  large  n.  To  investi- 
gate this  density  we  conducted  a simulation  of  44  points  uniformly 
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distributed  over  the  rectangle  with  sides  in  the  ratio  1;2,  which 
gave  rise  to  64  Delaunay  triangles.  These  numerical  values  corres- 
pond as  closely  as  possible  to  the  Iowa  data.  Thiosseii  polyqt>ns  woic' 
constructed  using  the  program  of  Green  and  Sibson  (1077).  we  will 

examine  this  hypothesis  in  §2.2.1  and  §2.2.2  through  the  dens i ty  fmclion  of 

A,  which  is  twice  the  area,  and  through  the  marginal  density  of  the 
minimal  angles. 

2.2.1.  Verification  through  Area 

For  reasons  given  later,  we  look  into  the  behaviour  of  the 
distribution  of  areas  of  Delaunay  triangles  of  unit  radius.  The 
areas  are  given  by 

JjLsin  20^+  sin  2a^-  sin  2 ( a^+  a^)  ] = *5A,  a^>  0,  a^>  0,  a^+  a^<  n (2.4) 

The  distribution  function  of  A , H(A),  from  f(a^,a2)  given  by 
(2.2),  is  found  to  be 

-1 

H(A)  = (3n)  [s(a)-  a)  sin  2a  da 

a=0 


where 


-2(3n) 


-1 


n 

J. 


sin  a sin  s(a)  da. 


a=0 


0<  A < 


3^ 

2 


(2.5) 


and 

s(a) 


cos  fmax(-l. 


sin  2a-A.  % 
2 sin  a'  ^ 


We  now  test  the  hypothesis, 

H : A ,...,A  are  independent  observations  from  H(A)  given  by 
o 1 n 

(2.5). 


I 
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There  will  be  some  form  of  dependence  on  the  A's,  but  it  has 
to  be  seen  whether  or  not  this  is  negligible.  Tabic  II  shows  a com- 
parison of  the  distribution  of  A for  the  siInul^^tcd  Poisson  process 
of  44  points  with  the  simulated  distribution  of  A under  Milos  den- 
sity, equation  (2.2). 

TAXU-K  II 

The  distribution  of  A(=2xAreii)  for  a simulated  Poisson  process  of 
44  points  KlvlnR  rise  to  64  Dolaunay  triangles  (obscrvell)  and  for 
the  simulated  values  from  the  Miles  distribution  given  by  (2.5), 
calculated  from  10,000  simulations 


Range 
of  A 

Observed 

frequency 

Expected 

frequency 

0.0 

-0.2 

1 

1.11') 

0.2 

-0.4 

1 

2.2.5 

6.50 

0.4 

-0.6 

0 

3.14  J 

0.6 

-0.8 

4 

3.67  ) 

8.04 

0.8 

-1.0 

8 

4.37  * 

1.0 

-1.2 

7 

5.10 

1.2 

-1.4 

7 

5.16 

1.4 

-1.6 

9 

5.38 

1.6 

-1.8 

8 

6.23 

1.8 

-2.0 

3 

6.30 

2.0 

-2.2 

2 

6.76 

2.2 

-2.4 

7 

7.00 

2.4 

-2.6 

7 

7.44 

Using  the  goodness  of  fit  criterion,  it  is  found  that  x^~  14.49. 

2 

The  5%  value  of  Xq  “ 16.92  and  hence  we  accept  the  null  hypothesis, 

Furtha;  it  is  found  that  for  the  simulated  Poisson  process  the  mean 

2 


and  the  variance  are  A = 1.53, 


S = .3665.  It  can  be  shown 
A 


( 

S 


r 
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frora  (2.2)  that 


, X tt  2 35  II 

^ = E(A)  = - , a = var(A)  = JJ  " "J" 


2.6] 


Under  the  null  hypothesis  for  large  n,  A ~ N(u,a  /n) . in  fact 
_1  k — 

a n lA-ul  = .48  and  consequently  we  again  accept  the  null 
hypothesis. 

2.2.2.  Verification  through  Marginal 
Consider  the  order  statistics  a, 

1 ; V ^ 

It  can  easily  be  shown 
that  the  marginal  density  of  minimum  angle,  becomes 


■(!)•  “(2)-  “-(3) 


“(1)''  “(2)''  “(3) 


“(D*  “(2)^  “(3)=  "■ 


f(a„  ) = (2/n)[(n  - sin  cos  2a(j)-  ““(11 


TJ 


° ^(1)^  3 • 


We  now  test  the  hypothesis, 

H : 0.,, are  independent  observations  from  (2.7). 
o (1)1  (l)n 

2 2 

Using  the  x goodness  of  fit  criterion,  x “ 8.18.  The  5%  point 


(2.7) 


2 

for  x^  ~ 14.07  and  so  again  we  accept  the  null  hypothesis. 


TABLE  III  (page  10) 

The  distribution  of  the  minimum  angle,  ® simulated  Poisson 

process  of  44  points  giving  rise  to  64  Delaunay  triangles  (observed) 
and  for  the  expected  values  for  Mile's  density  given  by  (2.7) 


. I 
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Kungo 


Observed 


of 

a 

(1) 

f requ' 

5° 

0 

- 

0 ( 
o 

1 

15 

3 

1.5*’ 

0 

o 

8 

o 

0 

25° 

13 

25° 

30° 

9 

50° 

35° 

9 

as"" 

o 

0 

6 

40° 

45° 

7 

45° 

- 

0 

o 

4 

c 

o 

LO 

55° 

4 

55° 

0 

o 

0 

Kxpeclod 

frequency 


.97 

2.88  I 

4.C>(> 

6.22 

7.13 


1.40 


10. 18 


For  the  simulated  Poisson  process  .53.  From  (2.7),  it 

can  be  shown  that 

, _ 27  2 , ^ _ V3  27  1 ,27 

^ ^ “"(I;  16n  ' var(a^^j)  n * 32  8 ~ ^16tt^  ’ 

Under  N{v,T^/n)  a.nd  so  we  have  t n ^ - v ] = , 22  . 

Hence  we  again  accept  the  null  hypothesis. 

Boots  (1974)  has  suggested  using  the  marginal  distribution  of 
a random  angle  from  a random  Delaunay  triangle  for  verifying  whether 
the  pattern  is  random.  He  uses  Miles  marginal  density 


f(a)  = 4f(77-a)  cosa  + sin  aj  , 0 < a < rr  , 


(2.9) 


and  suggests  randomly  selecting  a sub-sample  of  triangles  and  then 
choosing  one  angle  at  random  from  within  each  triangle.  This  pro- 
cedure would  of  course  throw  away  a substantial  part  of  the 


available  information. 


To  reach  a firmer  conclusion  concerning  the  null  hypothesis 
extensive  simulations  are  necessary,  however  on  the  basis  of  the 
one  simulation  conducted,  we  have  no  grounds  to  reject  H 

3.  VON  MISES  ALTERNATIVE 

For  any  triangle  there  is  a unique  circle,  the  circumcircle, 
which  can  be  drawn  to  pass  through  all  three  vertices  of  the  tri- 
angle. Using  this  fact,  we  can  view  the  vertices  of  a particular 
Delaunay  triangle  to  be  observations  on  the  circuradisc.  We  now  make 
the  following  assumptions: 

(I)  Central  Place  Theory  predicts  equilateral  Delaunay  triangles. 
We  propose  to  investigate  the  theory  by  considering  the  shape  rather 
than  the  size  of  the  triangles.  Considering  the  size  of  the  triangles 
to  be  unimportant,  each  triangle  is  scaled  so  that  it  is  circumscribed 
by  the  unit  circle,  i.e.  R = 1. 

(II)  With  the  angular  observations  thus  formed,  independence 
between  triangles  is  assumed.  We  have  already  given  some  justifica- 
tion for  this  assumption  in  §2. 

3.1.  The  von  Mises  Model 

Consider  ®i'®2'^3  Mises  variables  with  common 

concentration  parameter  K , then  9^  is  distributed  as  M(p^, K) 

for  i=l,2,3  and  therefore  the  joint  p.d.f.  of  9^,  9^,  0^  is 


(2n)^I^(it) 

o 


K[cos(9  - d.  ) +COS  (9  -p  )+cos{9  -M  ) ] , 
e 1 1 2 2 3 3 


tC  > 0 . 


(3.1) 


If  we  take  p*  = (0, 2n/3,4Ti/3)  it  is  obvious  that  as  K ♦ ® , the 
triangles  tend  to  become  equilateral  and  therefore  K can  be  regarded 
as  a measure  of  the  degree  of  equilateral ity. 
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3.2.  von  Mises  Arcs 


AS  the  model  stands,  it  cannot  be  directly  applied  to  Lhe  prob- 
lem as  we  do  not  know  the  individual  marginal  observations,  only  the 
angle  differences.  However  the  data  in  Table  I can  be  viewed  in  terms 


of  the  'arc  lengths'  2a^,  ^'^2’  ^*^3  unit  circle.  Therefore  we 


investigate  the  distribution  of  these  arcs  which  we  denote  by 


i , and  $ such  that  $ + “f  t = 2n,  and  we  can  write  'i’  , ^ in 

^2^  123  123 


terms  of  possibilities. 


(i)  If  ®i'®2'^3  anticlockwise  order 


$ = (9  - 9 )mod  2n,  $ = ^^-3“  6_)rood  2n,  $ = 9 ) mod  2n 


or 


(ii)  if  9,, 9 ,9^  are  in  clockwise  order 
12  3 


(9^-  92)mod  2tt,  ^ (6^-  8^)mod  2n,  (9^-  0^)mod  2rr. 


Consider  (3.1).  By  changing  variables  from  (9j^,92*9^)  to 

( (92“  9^)mod  2n,  6^,  ^ integrating  out  9^,  we 

arrive  at  the  joint  p.d.f.  (degenerate)  of  $ ^ . The  p.d.f.  of 

■L  f ^ 


($  .$  ) is  found  to  be 
1 2 


f(^l,'i'2)  = [(2n)^I^(K)]~^[l^fK[3+2cos($^-P2+t-‘j^)  + 2cos  ( $ ) 


+ 2cos  (^^-pj^+U^)  3 


+ I f k[3+2cos ( $ -p, +p_) + 2cos($_-P-+P-)+2cos($  -p  +p  3 }], 

O 112  22J  JJl 


iji  = 2n-$  -$  , 0<  $ s 2n-$  . 0<  2tt  . 

3 12  2 1 1 


(3.2) 


J 
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For  equilateral  triangles,  we  take  u’  = {0 , 2ti/3 , 4ti/3 ) and 
hence  the  model  becomes 

" [(2n)^I^(K)]"^[l^fK[3+2cos(#^-  ^)+2cos(^2-  ~) +2cos 

+ r {tc[3+2cos(‘l>^+^)+2cos($2+  — )+2cos('l'.^+  . 

= 2n-$^-$^,  0<  2n-$^,  0<  2ii  . (3.3) 

Particular  Cases 

For  K = 0,  the  density  collapses  to  become 

^ l/(2n^),  0<  2n-$^,  0<  2n  , (3.4) 

which  is  what  one  would  expect  from  the  construction.  For  small  K 
we  can  approximate  the  density  by  expanding  the  Bessel  functions  and 
ignoring  0(k‘^)  and  higher  powers  of  K which  gives 

2-1  ^ 

= (2n  ) [1-^  ) cos$.}  . 

^ ^ 4 j4^  3 

^2  = 2n-$^-$^,  0<  2n-$^,  0<  ^ 2Tr  . 

so  small  that 
Defining  6 = 

of  tc,  (3.3)  becomes 


For  large  K , $ . -*  ——  + 6.  where  6.  is 
1 3 1 1 

3 ^ 

0(6.).  We  have  therefore  > 6.=  0. 

^ i4l  " 


neglecting  o(K  ^)  and  smaller  powers 


(3.5) 

we  can  neglect 
7^6^  and 

^ i4l  " 


1 


I 


( 

I*. 


I 
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2n^/3 


-J5K  6 


::2Ji<6  ^ IIL  _ 6 
3 2 3 1 


-21T 


6,  -i  — 


4_^ 

3 


Hence  ^ bivariate  normal  distribution  with 

2 2 

^2“  (2n)/3,  a^=  a^=  2/k  and  p = -1/2.  This  same  result  could 

have  been  found  by  taking  the  normal  approximation  to  the  von  Miscs 

density  in  (3.1)  and  using  the  fact  that  the  difference  of  two 

normal  variables  is  again  normal. 

We  now  investigate  the  behaviour  of  the  distribution  of  A and 

from  this  model  for  varying  K.  Since  the  exact  form  of  the 

densities  is  unknown,  we  turned  to  simulations  of  the  model.  The 

method  of  simulation  employed  was  that  of  the  Rejection  Procedure. 

For  each  value  of  k,  ®ii'®2i'^3i  simulated  independently  and 

the  corresponding  values  of  A.  and  . were  calculated  for 

1 (1)1 

i=l, . . . , 10, 000.  The  results  are  presented  in  Tables  IV  and  V (p.l5). 

Table  IV  shows  the  density  of  A changing  from  an  almost  mono- 
tonic decreasing  function  at  (C=0,  to  an  almost  monotonic  increasing 
function  at  k=2.8.  increasing  K still  further,  results  in  the 
density  becoming  heavily  concentrated  around  equilaterality  and 
being  almost  entirely  contained  in  the  upper  half  of  the  region. 

In  Table  V a similar  process  is  occurring,  but  here  the  mode  is  not 
given  by  the  maximum  value  of  ^(i)* 


Identification  with  Miles'  Density 

The  question  arises  whether  the  von  Mises  density  is  related  to 
Miles'  density  for  some  K.  To  this  end,  consider  the  area  of  a 
particular  triangle  which  is  given  by 


I 

( 


ilau);)! 
of  A 

K = 0 

K = 2.8 

. ) 

< - 12  i 

0.0 

-0.2 

2 124 

.12.5 

O 

‘>.2 

- 0.4 

1150 

3 1 8 

O 

- 0.(; 

802 

:i4M 

o 

0 .G 

- O.H 

802 

372 

'I 

0.8 

- 1.0 

757 

4G4 

^ 1 

1 .0 

- 1.2 

GOO 

547 

8 

1 .2 

- 1.4 

552 

G12 

10  • 

1 

1 .-1 

- l.G 

503 

G89 

43  ■! 

1 .6 

-1.8 

525 

800 

140  j 

1.8 

- 2.0 

525 

1025 

382  ' 

2.0 

- 2.2 

502 

1283 

961  !] 

2.2 

- 2.4 

471 

1480 

2427 

2.4 

-2.0 

408 

1647 

6012  |i 
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Area  = *2A  = b 1 sin  ( 0 ^ -0 ^ ) +s  in  ( 0 ^-6^ ) +sin  ( 6^ -B ^ ) 1 , 

Area  has  a certain  intuitive  appeal  in  that  for  a given  cir- 
cumcircle,  the  triangle  with  the  maximum  area  is  an  oqui lateral 
triangle.  The  modulus  sign  makes  A awkward  for  an.)  I yt.iccii  pur- 
poses, so  we  consider  instead  A^  which  has  the  same  property  of 
being  a maximum  for  equilateral  triangles.  Hence  we  identify 

E(A  ) for  both  densities  and  it  is  found  that  K , the  value  of 

o 

K corresponding  to  Miles'  density,  should  satisfy 


^5 

12 


3 3 2 


(3.7) 


where 


A (K  ) 
r o 


The  solution  of  (3.7)  gives  K = 2.20. 

o 

We  now  investigate  how  close  this  approximation  is  in  fact. 

This  is  done  through  the  distribution  of  A and 

Miles'  density  and  the  von  Mises  model.  Since  the  exact  forms  of 
the  distributions  are  not  known,  except  for  under  Miles' 

density,  (2.7),  most  of  the  comparison  is  done  through  similation. 
AS  before  the  method  employed  is  the  Rejection  Procedure.  The 
results  are  shown  in  Figs.  2 and  3. 


As  the  Figures  show,  the  two  distributions  agree  fairly  well 
for  both  A and  except  in  the  lowest  region  whore  the  von 

Mises  model  predicts  a much  greater  number  of  small  triangles. 

This  particular  point  is  discussed  in  the  next  section. 

3.4.  Appraisal  of  the  von  Mises  Model 

Under  the  usual  circumstances,  the  test  for  randomness  with  the 

von  Mises  density  is  H : K = 0 V h.  : K > 0,  however  this  is  not 

o 1 

applicable  here.  If  the  asymptotic  approximations  of  Miles  are 

valid  for  n as  small  as  44,  which  the  simulation  described  in 

Table  II  would  appear  to  confirm,  the  test  becomes  of  the  form 

H ;K=K  VH,  :K>K  where  K = 2.20. 
o o 1 o o 

It  has  to  be  strongly  emphasized  that  this  approach  is  not  de- 
scribed through  a spatial  process  and  in  that  sense,  it  has  to  be 
admitted  that  the  procedure  is  'ad  hoc'.  For  large  K the  model 
has  certain  good  features  (e.g.  sensible  m.l.e.'s)  and  it  is  cer- 
tainly a viable  working  alternative,  especially,  once  having  re- 
jected the  hypothesis  of  randomness.  For  K = it  can  be  thought 
of  as  similar  to  Miles,  but  for  small  K,  K < K^,  the  density  appears 
to  be  meaningless  in  the  spatial  content  and  so  can  be  said  to  be 
non-spatial.  This  occurs  due  to  the  relatively  large  probability 
of  colinearity  under  the  von  Mises  model. 

Under  the  Miles  density,  equation  (2.2),  f(a^,a2)  disappears 

when  the  angles  are  colinear,  whereas  since  I (0)  = 1,  (3.3)  is 

o 

still  equal  to  a positive  constant,  and  so  is  bounded  away  from 
zero.  As  K increases  this  constant  tends  to  zero,  but  for  k<k 

o 

it  is  quite  significant. 

In  an  actual  realisation  of  a spatial  process  for  a moderate 
number  of  points,  for  the  three  points  of  a Delaunay  triangle  to 
be  almost  colinear,  the  circumdisc  of  the  triangle  would  be  very 
large.  The  alternative  definition  of  a Delaunay  triangle  is  that 


3 points  form  a Delaunay  triangle  iff  their  circumdisc  is  empty. 
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and  so  given  3 points  practically  colinear,  by  the  very  definition 
of  a Delaunay  triangle,  it  is  most  unlikely  that  the  3 points  would 
form  such  a triangle. 

It  would  be  nice  if  we  could  arrive  at  a von  Mises  type  density 
through  a non-Poisson  process,  but  this  does  seem  to  be  out  of  out- 
reach at  the  present  time,  if  it  is  indeed  at  all  possibU'. 

4.  ESTIMATION  FOR  THE  VON  MISES -TYPE  MODEL 
It  is  reasonable  to  assume  that  the  p's  are  known  and  for 
which  we  take  p'=  (0, 2n/3,4TT/3)  . Consequently  we  need  to  estimate 
the  concentration  parameter  IC.  Two  methods  of  estimation  are 
given. 

4.1.  Method  of  Moments 

Defining  0^,  independent  von  Mises  variables  with 

p'=  (0, 2n/3,4n/3)  and  common  concentration  parameter  k:,  it  is 

found,  as  in  (3.7)  that 


E(A^)  = f + + 3aJ{K)A2(IC)  +fAj(K)  , 


(4.1) 


where  A^(K)  is  defined  in  §3.  We  take  as  our  estimate  of  E(A  ), 


A 2 1 

E(A  ) = - 


- y ^ 

" ik  ^ 


(4.2) 


where  Aj^,...,A^  are  the  observed  areas  of  the  triangles.  On 
equating  (4.1)  and  (4.2)  we  can  obtain  the  moment  estimate  K of  tc 
Unfortunately  the  equation  can  only  be  solved  numerically. 

Example 

From  Table  I we  can  calculate  (4.2)  for  the  Iowa  data,  which 


gives 


V 

n ) A. 


. = 3.4073.  Substituting  into  (4.2)  and  solving 


numerically,  it  is  found  that  K = 2.80. 
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4.2.  Maximum  Likelihood  Estintation 

Let  ^2i*  ^3i^'  i-l»...»n  be  the  observed  arc  longtlis. 

From  (3.8)  it  can  be  easily  shown  that  the  maximum  likelihood 
estimate  of  K,  is  given  by 


A 

I^(K) 

~A~ 

lo(^) 


(4.3) 


, A A 

I (Ku. ) + I (Kv. ) 
or  or 


where  u = [3+2cos($  -•^)  + 2cos($  -^)  + 2cos($  .-^)]^ 

1 Xx^  2X-J  3i3 


and  v^  = [3+2COS  ( 4-^^+ -^  ) + 2cos(4>^^+^)  + 2cos  ( $^  . + ^) 


By  the  methods  of  §3.2  we  can  find  approximations  to  the 

maximum  likelihood  equation  for  small  and  large  values  of  K. 

A 

For  small  K,  by  expanding  the  Bessel  functions  and  neglecting 
A4 

0 ( It  ),  (4.3)  becomes 


(4.4) 


where  C . = ) cos  ♦ . 

1 j4l  31 


For  large  K,  define 


< * .ill  - 1 

6..=  . 6=r- 

11  j 1 3 3n 


1^2/37  i 


Following  §3.2,  (4.3)  then  reduces  to  K = 2/6.  But  also  we  have 
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n 


A = - y A. 


„ 3V3 


(2  - 6).  Hence, 


K = (3V3)/(3V3  - 2A)  . 


(4.5) 


A 


Consequently,  K depends  only  on  a for  large  K. 
Example 


Consider  i~l»****63»*  i=l,2,3  where  the  a's 


are  given  in  Table  I for  the  lowa  data.  Substituting  directly 


into  (4.3)  and  solving  numerically  gives  K = 2.83,  which  is  in 
good  agreement  with  our  estimate  of  K,  K = 2.80,  given  by  the 
method  of  moments. 

A 

Now  let  us  consider  the  large  K approximation.  it  is  found 
that  A = 1.7105.  Substituting  this  value  of  A into  (4.5) 


gives  the  m.l.e.  of  K as  2.93,  which  is  only  slightly  larger 


than  K = 2.83. 


5.  ANALYSIS  OF  THE  IOWA  DATA 


We  can  now  analyze  the  Iowa  data  given  in  Table  I. 
5.1.  Models 


We  first  enquire  whether  the  data  could  be  described  by  a 
von  Mises  type  model  with  p'=  (0, 2n/3,4n/3) . This  we  do  through 


the  marginal  distribution  of  A and  We  note  from  §4.2 

A 


that  the  m.l.e.  of  K is  K = 2.8.  Tables  VI  and  VII  give  a 
comparison  between  the  Iowa  data  and  the  expected  frequencies  for 


the  von  Mises  K =2.8  for  A and  *(i)'  given  by  10,000  simu- 


lations . 


J 
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, TADI.K  VI 

I Comparison  of  A for  Iho  observed  Iowa  data  and  Mioso  values 

I expected  from  10,000  simulations  for  tlui  von  Mlsos  model  wltb 

K = 2.8  for  Miles'  density 


Hatige 

Observed 

E 

E 

2 

of 

A 

frequency 

1 

0.0  - 

O.G 

1 

C .21 

6.42 

0.6  - 

1 .0 

7 

5.26 

7 .93 

1.0  - 

1.4 

14 

7.31 

10.11 

1.4  - 

1.8 

7 

9.95 

11.44 

1.8  - 

2.2 

13 

14.54 

12.87 

2.2  - 

2.6 

21 

19.70 

14.23 

■ 


I 

i 

I TAllbK  VII 

1 Compurlnon  of  values  between  the  observed  Iowa  data,  those 

’ (1) 

^ values  expected  fur  the  von  Mines  model  with  k = 2.H  given 

j by  10,000  simulations  and  tiioso  values  expected  for  Miles' 

j ' density  (E,,) 


Hungc  Observed 

of  frequency 


0® 

- 20° 

0 

5.35 

3.80 

20** 

- 40° 

8 

7.37 

10.71 

“^O 

60 

- 60° 

14 

11  .65 

15.37 

- 80° 

16 

15.90 

16.17 

HO^ 

loo” 

- 100° 

18 

15.99 

12.32  1 

- 170° 

7 

6 .71 

4 .62  i 

i 
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using  a X goodness  of  fit  test  for  A and 


testing  against  ~ 9.49  and  X4  < ~ 13.28  in  both 

2 

cases,  it  was  found  for  A,  X ~ 12.22  and  for  i 


(1)  ' 


0.13. 


Hence  we  just  reject  the  fit  for  A at  the  5%  level,  but  accept 

the  fit  of  the  model  for  ^(i)* 

The  final  columns  in  Tables  VI  and  VII  give  the  corresponding 

frequencies  under  Miles'  model  from  simulation.  Applying  the 

test  for  Miles'  density  produces  the  mildly  surprising  result  tliat 

we  would  accept  the  fit  of  this  distribution  to  the  Iowa  data. 

2 2 

For  A we  have  x ~ 10.46  against  X5  ( 5%)  = 11.07  and  for 

2 2 ' 

X = 6.87  against  X3(^%)  “ 7.81.  The  best  procedure  would  have 

been  to  look  into  the  bivariate  distribution  of 

sample  size  of  n = 63  is  anyhow  too  small  to  draw  clear-cut 

conclusions.  However,  as  we  have  seen  that  for  K =2.2  the 

o 

von  Mises  type  model  is  adequate.  We  will  study  the  problem  more 
fully  from  this  angle  in  the  next  section. 


5.2.  central  Place  Theory 

We  wish  to  investigate  whether  the  Delaunay  triangles  of  the 

Iowa  data  can  be  regarded  as  having  arisen  as  the  result  of  a 

Poisson  process  in  the  plane  or  whether  they  occurred  through  the 

action  of  a more  regular  force.  Under  the  von  Mises-type  model 

of  §3,  this  problem  can  be  looked  upon  as  testing  H^  : K = 

against  H,  : IC  > K where  p'  = (0, 2n/3,4n/3)  . if  X is  the 
1 o — 

likelihood  ratio  for  the  problem,  it  is  found  that 

n 

-2  log  X = 6n  log  (K  /K  ) +27  log  ([l  (K  u.)  + I (K  v )]/ 

O X . 0X1  0X1 

1^ 

[l  (K  u.  ) + I (K  V.  ) ]) 

001  001 
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where  tc  = 2.20,  K = 2.83  and  u.  and  v.  are  given  by  (4.3). 
Substituting  the  values  from  Table  I gives  -2  log  X - 4.08. 

2 

Xj^(5%)  = 3.84  and  hence  we  reject  K = at  the  5%  level. 

With  the  evidence  before  us,  there  seems  justifiable  grounds  at 
the  5%  level  to  suppose  that  the  Iowa  data  is  a result  of  a pro- 
cess that  is  more  regular  than  the  Poisson  and  hence  there  is 
some  evidence  that  central  Place  Theory  holds  for  this  particular 
region. 
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ABSTRACT 

Central  Place  Theory  predicts  a regular  spatial  pattern  in 


the  plane  and  we  observe  that  the  Delaunay  triangles  will  be 
equilateral  under  the  theory.  However,  when  the  pattern  is 
'random'.  Miles  (1970)  has  given  the  asymptotic  p.d.f.  of  the 
interior  angles  of  a random  Delaunay  triangle.  We  propose  a 
von  Mises-type  model  with  a concentration  parameter  K;  the 
larger  the  value  of  K , the  closer  we  are  to  the  Central  Place 

y 
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Theory.  We  show  that  the  model  can  be  approximated  to  th<^  Miles' 
density  for  some  value  of  K.  we  provide  the  moment  and  maximum 
likelihood  estimators  of  K,  and  it  is  recognized  that  the  areas 
of  the  Delaunay  triangles  play  an  important  role.  We  construct 
a test  of  departure  from  the  random  pattern  with  the  alternative 
of  Central  Place  Theory.  As  a numerical  example,  we  analyze  44 
Central  Places  in  Iowa  where  we  find  some  evidence  for  the  valid- 
ity of  Central  Place  Theory  in  that  particular  region. 

RESUME 

La  Theorie  de  la  Position  Centrale  envisage  une  repartition 
reguliere  des  implantations  urbaines  dans  la  plaine.  Dans  cette 
theorie  on  s ' attend  a ce  que  les  triangles  de  Delaunay  de  la 
repartition  soient  equilateraux.  Cependant,  quand  la  repartition 
est  "au  hasard".  Miles  (1970)  a donne  la  distribution  asymptotique 
pour  les  angles  interieurs  d ' un  triangle  de  Delaunay  pris  au 
hasard.  On  nous  propose  un  modele  de  von  Mises  avec  un  parametre 
de  concentration  "K".  La  validite  de  la  theorie  de  la  position 
centrale  augmente  avec  "K".  on  demontre  que  1 ' on  peut  donner  une 
approximation  du  modele  par  la  densite  de  Miles  pour  uxe  valeur 
donnee  de  "K".  on  trouve  des  evaluations  de  "K"  en  utilisant  les 
moments  et  1' equation  du  maximum  de  vraisemblance . L' importance 
des  surfaces  des  triangles  de  Delaunay  est  reconnue,  et  une 

X * 

experience  entre  la  "priseau  hasard"  et  la  theorie  de  la  position 
centrale  est  au  cours  de  developement . A titre  d'exemple 

* s 

numerique,  on  analyse  44  positions  centrales  dans  1 ' Iowa  ou  1 ' on 
trouve  des  preuves  de  la  validite  de  la  theorie  dans  cette  region 
particuliere . 
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